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Abstract: The CDMS Collaboration has reported two candidate events for dark matter. 
If the events are due to the elastic scattering of dark matter, the dark matter would be a 
WIMP dark matter with its mass of the order of 10-100 GeV and its scattering cross section 
with a nucleon is about 10~'^^-10~''^ cm^. We show that such a dark matter is properly 
realized as a neutralino dark matter in the light higgs boson scenario of the MSSM. The 
lightest higgs boson mass can be lighter than 114.4 GeV in the scenario because of a 
suppressed interaction between higgs boson and Z bosons. As a result, a large scattering 
cross section between the dark matter and ordinary matter is obtained. 
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1. Introduction 

Existence of dark matter in our universe has aheady been established thanks to recent 
cosmological observations such as the Wilkinson Microwave Anisotropy Probe (WMAP) 
experiment |jl| , and now we know that about 23% of the energy density of the universe is 
provided by dark matter. We know, however, very little of what the dark matter is. Many 
attempts to identify the dark matter have been performed so far [Q]. Very recently, two 
candidate events for dark matter have been reported by the Cryogenic Dark Matter Search 
(CDMS) collaboration The signals may be due to the scattering between a nucleus in 
the detector and dark matter in the halo associated with our galaxy. Though the events 
are still consistent with background fluctuation at the probability of 23%, the events may 
give us important information about the dark matter. 

If the two events are caused by the elastic scattering of the dark matter inside the 
detector, the dark matter would be a Weakly Interacting Massive Particle (WIMP) with 
its mass of 10-100 GeV [^. It is also important to notice that the scattering cross section 
between dark matter and a nucleon is of the order of 10~^^-10~^^ cm^ if the scattering 
occurs through spin-independent interactions. Such a large cross section can be obtained 
from the process in which the higgs boson is exchanged, when the interaction between 
the higgs boson and two dark matters is large enough. This kind of scenario has already 
been discussed by many papers even in a very short period after the report of the CDMS 
collaboration |5|. 

In this letter, we propose another scenario to account for the CDMS result, where 
the large scattering cross section is obtained by a light higgs boson with mass less than 
114.4 GeV. One might think that such a light higgs boson is not favored by the LEP2 
experiments |^] . However the higgs boson indeed can be light in the framework of the light 
higgs boson scenario (LHS) of the minimal supersymmetric standard model (MSSM) 0, 
|8|, |9|, |l^. Furthermore, as shown in Ref. Q, the scenario is consistent with the WMAP 
experiment as well as other experimental constraints such as 6 ^ 57 and Bg — )• 
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The essential reason why the higgs boson can be hghter than the LEP bound {rrih > 
114.4 GeV) in the LHS is as follows. There are two higgs doublet fields in the MSSM, which 
leads to two CP-even neutral higgs bosons. By choosing an appropriate mixing between 
these two, the couphng constant between the lightest higgs boson and Z bosons can be 
significantly smaller than that of the standard model (SM). The mass of the lightest higgs 
boson therefore can be lighter than the LEP bound, because the bound comes from the 
process e^e~ — ?> Zh. In addition, the LHS can explain the 2.3 a level excess of the events 



at 98 GeV in the LEP2 experiments ^, IC], which is very difficult to be explained in the 
SM. 

We consider the SUSY model with non-universal scalar masses for the higgs multiplets 
(NUHM), which is one of concrete models realizing the LHS being compatible with the 
grand unified theory (GUT). Details of our setup will be shown in the next section. In order 
to clarify the region consistent with many experimental constraints in a broad parameter 
space, we use the Markov Chain Monte Carlo (MCMC) method, which will be discussed in 
section ^. Finally, in section it is found that the model predicts the bino-like neutralino 
dark matter with its mass 50-250 GeV. Interestingly, the scattering cross section between 
the dark matter and a nucleon can be around 10"^'* cm'^, which is nothing but the value 
requested by the CDMS result. 

2. LHS of the MSSM 

With the assumption of CP conservation, two higgs doublet fields {Hu,Hd) in the MSSM 
leads to the spectrum of the higgs sector composed of three neutral and two charged scalar 
bosons: two CP-even higgs bosons (h and H, where h is lighter than H), CP-odd higgs 
boson (A), and charged higgs bosons (H^). The mass eigenstates of the CP-even higgs 
bosons are given by the mixing states between neutral components of and -ff^. 



sin a cos a\ /Re 



cos a sin a \ Re H. 







(2.1) 



where a is the angle introduced to diagonalize the mass squared matrix of CP-even higgs 
bosons, 



m\sl + mlc^p + Add ~{m\ + m%)spCj3 + Adu 
-{m\ + m\)spci3 + Adu m^c^ + m|s^ + A„„ 



(2.2) 



The mass of the CP-odd higgs boson {Z boson) is denoted by tua {mz) and Cji{sp) = 
cos /3 (sin /3) , where the ratio of the vacuum expectation value of the higgs doublet fields is 
given by tan/3 = {H^)/{H^). The radiative correction to each component in the matrix is 



represented by Ajj, where its detailed expression can be found in Ref. [11] 



The higgs boson has been searched using the process e^e~ — ?• Zh in the LEP2 experi- 
ments. The coupling constant between h and Z bosons is given by gzzh = S'^Zh sin(/3 — a) 
in the MSSM, where g^^^l is the corresponding coupling constant in the SM. As a result, 
the lightest higgs boson h can be lighter than the LEP bound (m/j > 114.4 GeV) when the 
coupling constant gzzh is significantly smaller than g^^^jl 0, |8|, ||, . 
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Since the angle f3 is expected to be as large as 7r/4-7r/2, a large mixing angle a is 
required to realize a small gzzh- Such a large angle is obtained when tjia is small enough, 
which leads to the situation that all higgs bosons in the MSSM are at the scale of O(IOO) 
GeV. Moreover, suppressed sin(/3 — a) leads to enhanced cos(/3 — a), which results in a large 
coupling constant gzzH between heavy higgs boson and Z bosons. Therefore, the heavy 
higgs boson should satisfy the constraints of the LEP higgs boson search. In addition, 
the coupling constant gzAh between the lightest higgs boson, CP-odd higgs boson, and Z 
boson is also proportional to cos(/3 — a). This coupling leads the e+e" — )■ hA process at 
the LEP2 experiments. Hence, the constraint for this process should also be taken into 
account. Fortunately, this constraint is not severe due to the P-wave suppression, because 
the coupling gzAh originates in a derivative interaction. 

The LHS cannot be realized in the constrained MSSM, which is widely used to study the 
MSSM. If the universal masses of sfermions and gauginos are fixed to derive < 114.4 
GeV, the masses of other SUSY particles becomes too small. In order to realize the 
LHS, the masses of the higgs doublet fields should be different from others as ^^^^ = 
(1 + Shu{Hci))^o- It is reasonable because the higgs multiplets are not necessarily in the 
same multiplet of GUT with other scalar particles. The simplest model with this boundary 
condition is the NUHM. Hence, we adopt this model as a reference model to investigate 
the LHS iQ. The NUHM has six free parameters, (mo, ?7ii/2! ^O; tan/S, /i, m^). The first 
three parameters (mo, mi/2, ^o) are universal sfermion mass, gaugino mass, and tri- linear 
coupling, which are defined at the GUT scale Mq- Other parameters, tan/3, higgsino mass 
//, and rriA are defined at the electroweak scale. This parameterization allows us to treat 
the masses of two higgs doublet fields at Mq as free parameters. Using these values at the 
mz scale, a boundary condition at Mq is derived by the renormalization group running. 



Masses of SUSY particles then run back from Mq to mz- We use ISAJET 7.75 [13| to 
evaluate the renormalization group running. 

3. Markov Chain Monte Carlo 



Markov Chain Monte Carlo method is a random sampling algorithm, which constructs 
a series of parameter sets (Markov chain) |jl^. The samples of the chain obey a distribu- 
tion which is proportional to a given distribution function. The distribution function of 
our interest is a posterior probability distribution function of model parameters x under 
experimental data D. Bayes' theorem tells us that the posterior probability distribution 
function P{x\D) satisfies the following equation, 

P(x\D) - ^(^1^)^^^) (3 1) 

where P{x) is the prior probability function reflecting our knowledge about the model 
parameters x, while P{D\x) is representing the likelihood for the distribution function of 
experimental data D at given model parameters x. In our analysis, a linearly flat prior has 
been used for P{x), where P{x)dx gives a constant probability. 
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Constraints 



References 



Relic abundance of dark matter {Q,DMh ) 

BR{b 57) 

/9-parameter (A/o) 

BR{Bs fi+fi-) 

Upper bound on gzzh {gzzu) 

Upper bound on gzAh {qzah) 

Lightest neutralino mass 

Chargino mass 

Right-handed selectron mass 

Right-handed smuon mass 

Right-handed stau mass 

Sneutrino masses 

Stop and sbottom masses 

Gluino mass 

Squark masses (1st and 2nd generations) 



0.1099 ±0.0062 

(3.52 ±0.25) X 10"^ 

± 0.0009 

< 5.8 X 10"*^ 

Function of {mn) 

Function of mh{mH) + rriA 

> 50.3 GeV 

> 103.5 GeV 

> 99.9 GeV 

> 94.9 GeV 

> 86.6 GeV 

> 94 GeV 

> 95 GeV 

> 308 GeV 

> 379 GeV 



(WMAP only) 



(p. 137) 
IE 






m 



m 



m 



m 



m 



m 



m 



m 



22 



Table 1: Constraints used in our Markov Chain Monte Carlo analysis. 



For experimental data D, we have used constraints shown in Table ||. The branching 
ratio BR(6 57) has been calculated with SusyBSG 1.3.1 [15]. Other observables and 



parameters are given by micrOMEGAs 2.2.CPC [16|. We then use a Gaussian distribu- 
tion for constraints on Qdm^^, BR(6 — )■ 57), and Ap, while other constraints have been 
introduced simply as boundaries of the model parameter space. In addition, since we are 
interested in the LHS, we concentrate on our exploration in the range m/i < 114.4 GeV. 
Finally, we have generated a chain of about 2.3 x 10'' samples. In figure |l|, the distribution 
of the samples is shown on the (m/j, sina) plane. Two distinct regions can be seen in the 
figure. In the region with low the property of h is essentially the same as the one in 
previous studies for rrih > 114 GeV. On the other hand, the high |a| region is nothing but 
the region discussed in section ^, which is characteristic of the LHS. We thus use only the 
samples in the high \a\ region by imposing a condition nih < 110 GeV in the following 
analysis, which is about 2.0 x 10'' samples. 



4. Direct Detection of Dark Matter 

Since all higgs bosons in the LHS are at the electroweak scale, detection rates at direct 
detection experiments of dark matter, which is the lightest neutralino in our analysis, are 
expected to be large ||^, ^, The distribution of the MCMC samples is shown in figure ^ 
with cells which have the same size in the linear scale. The right panel is the result with 
samples which satisfy also a condition Aa^ = (29.2 ± 2 x 8.6) x 10^-'^'^ @ (p. 482) for the 
anomalous magnetic moment [g — 2) of muon, while the condition is not used in the left 
one. The horizontal axis of these figures is the mass of the dark matter, and the vertical 
one is the spin-independent cross section between the dark matter and a nucleon. In our 
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Figure 1: The distribution of samples. The yellow (bright) region is of a large number of samples. 
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Figure 2: The prediction for direct detection experiments of the dark matter in the LHS. The 
figures show the distribution of samples which satisfy rrih < 110 GeV. The right figure is given by 
samples which satisfy an additional condition on g — 2 of muon. The yellow (bright) region is of a 
large number of samples. 



calculation, the y parameter, which characterizes the strangeness component in a nucleon, 
is set to be zero [0^]. Current bounds on the scattering cross section from XENONIO p6||. 



XENONIOO 1 27], and CDMS II 1^, ^ experiments are also shown in these figures as 
dashed, dot-dashed, and solid lines, respectively. Note that regions of a large number of 
samples in figure ^ do not mean the regions of good agreement with experimental results 
because the number of constraints (i.e. r^oM^^; BR(6 — )• sj), and Ap) is less than the 
number (six) of parameters in the NUHM. The cell with many samples simply shows that 
the values of the dark matter mass and the scattering cross section in the cell are easily 
obtained by the MCMC method because these values are predicted in a wide area of the 
parameter space of the NUHM. In figure ^, we also plot likelihood maximum^ which is 



""^We define the likelihood so that it is proportional to the product of Gaussian distribution functions for 



- 5 - 




^1 ^1 




Dark matter mass (GeV) Dark matter mass (GeV) 

Figure 3: Likelihood maximum which is obtained by calculating the maximum in each cell with 
samples which satisfy m/j < 110 GeV. The right figure is given by samples which satisfy an ad- 
ditional condition on g — 2 of muon. The yellow (bright) region is of a large value of likelihood 
maximum in a cell, where constraints on flj^Mh^, BR(6 57), and Ap can be satisfied in a good 
accuracy. 



obtained by calculating the maximum in each cell. The right panel in figure ^ is obtained 
by samples which satisfy a condition for g' — 2 of muon. A large value of the likelihood 
means that constraints in Table |^ are well satisfied. It is seen that most of the region of 
the samples agrees with the constraints in a good accuracy. The cross section is predicted 
to be around 10^^^-10~^^ cm^ in most of samples with the mass of the dark matter being 
around 100 GeV. If the CDMS II events are really coming from the elastic scattering of 
dark matter, the LHS can easily explain those events. By comparing the left and right 
panels in figures ^ and ^, it can be seen that a large part of the region is consistent with 
the result of — 2 of muon. 

The lower bound on the mass of the dark matter in figures |^ and ^ are controlled by 
the bound on the lightest neutralino in Table [l|, which is originally determined by those of 
charginos through the GUT relation on gaugino masses. The relic abundance of the dark 
matter around the lower bound is governed by the s-channel diagram to bb, in which the 
CP-odd higgs boson mediates. Such a region vanishes in right panels of figures |and|. On 
the other hand, the upper bound on the mass of the dark matter comes again from those 
of charginos. Too heavy charginos are not favored because the cancellation between the 
contributions from the charginos and light charged higgs bosons in the LHS is required to be 
consistent with the constraint on the 6 — > 57 process. In the region where the dark matter 
mass is larger than about 100 GeV, the relic abundance of dark matter is mostly governed 
by the processes whose final states are composed of higgs bosons |29| or the coannilation 



process between the LSP and stau; for samples of large scattering cross section around 
10""^^ cm^, the former processes are dominant, while small scattering cross section below 
around 10"'^^ cm^ is realized by the coannihilation process. The coannihilation process is 



i^umh^, BR(6 — i- sj), and Ap. The set of the central values for them in Table |^ makes the likelihood unity. 
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efficient in the samples used for tlie riglit panels of figures I and |. 

Even if the two events reported by the CDMS collaboration are due to background 
fluctuation, the neutralino dark matter in the LHS will be tested in near future at Super- 
CDMS, XENONIOO, and XMASS experiments Q because the region predicted by the 
scenario is not far away from the current bounds. 
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